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ABSTRACT 

We consider the precession rates of eccentric discs in close binaries, and compare the- 
oretical predictions with the results of numerical disc simulations and with observed 
superhump periods. A simple dynamical model for precession is found to be inade- 
quate. For mass ratios £ 1/4 a linear dynamical model does provide an upper limit for 
disc precession rates. Theory suggests that pressure forces have a significant retrograde 
impact upon the precession rate (Lubow 1992). We find that the disc precession rates 
for three systems with accurately known mass ratios are significantly slower than pre- 
dicted by the dynamical theory, and we attribute the difference to pressure forces. By 
assuming that pressure forces of similar magnitude occur in all superhumping systems, 
we obtain an improved fit to superhump observations. 
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1 INTRODUCTION 

Superhumps are now commonly found in short period cat- 
aclysmic variables (Patterson 1998) and X-ray binaries 
(O'Donoghue & Charles 1996). The most likely explanation 
is that these periodic luminosity variations are caused by 
the tidal stresses on an eccentric, precessing accretion disc 
(Whitehurst 1988; Hirose & Osaki 1990; Lubow 1991). In 
this model, the superhump period, Psh, equates to the period 
of the disc's precession, Pd, as measured in the binary frame. 
Although reliable measurements of Psh have been made for 
at least 53 separate systems (Patterson 1998), the observa- 
tions have yet to be compared properly with the model's 
theoretical predictions. 

In this paper we consider the factors determining the 
precession rate of an eccentric disc, and compare our the- 
oretical understanding, simulation results and the observa- 
tions. 



2 DISC PRECESSION PERIODS 

It is now well established that in binaries with mass ratios 
g ^ 1/4 it is possible for eccentricity to be excited in the 
accretion disc at the 3 : 1 eccentric inner Lindblad resonance 
(Lubow 1991). Here we define q — ra2/m\ to be the mass 
of the donor star divided by the mass of the accreting star. 
Under certain circumstances, eccentricity may be excited in 



systems with mass ratios as large as 1/3 (Murray, Warner 
& Wickramasinghe 1999). 

How an eccentric disc precesses coherently has been a 
source of confusion. For a single particle, the rate at which 
an elliptical orbit precesses 
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Now, for a given gravitational potential, a particle's angular 
frequency, f2, and radial or epicyclic frequency, k, are both 
functions of r. Hence the precession rate, uj, also depends 
upon the mean radius of the orbit. Yet neither the obser- 
vations (Patterson 1998) nor the simulations (Murray 1998) 
yield any indication of differential precession. So how does 
a disc, that extends over a range of radii, organize itself to 
precess with a unique frequency? 

The answer is that we are not dealing with a collection 
of isolated test particles, but with a gaseous disc. The eccen- 
tricity is excited at the Lindblad resonance, and then prop- 
agates inwards through the disc as a wave, getting wrapped 
into a spiral by the differential rotation of the gas (Lubow 
1992). If the spiral is wound up on a length scale much 
smaller than the radius (the tight winding limit), then u, 
the rate of azimuthal advance of the eccentricity, is governed 
by the dispersion relation 
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where c is the sound speed of the gas, and k is the radial 
wavenumber of the spiral (see e.g. Binney & Tremaine 1987). 
Now, as Lu is much less than fl and k we have 
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u) is determined for the wave in the region in which it is 
launched. As the wave propagates inwards, the frequency 
remains constant but the wavelength changes in response to 
the changing environment through which it moves. 

The above arguments are based upon the analysis of 
Lubow (1992). He also identified a third factor that con- 
tributed to the precession when the magnitude of the eccen- 
tricity was changing secularly. In this paper we are interested 
in the long term mean value for lu and so this extra term 
need not be considered further. 

Equation ^ tells us that a gaseous disc will precess 
more slowly than a ballistic particle at the resonance radius. 
Lubow (1992) estimated that pressure effects could reduce 
uj by approximately one percent of fi. Now ui^y^ is itself of 
the order of a few per cent oi Q., so the reduction is sig- 
nificant. (Simpson & Wood (1998) misinterpreted Lubow's 
work and ignored the pressure term as being a few percent 
of the dynamical precession). 

Hirose & Osaki (1993) obtained estimates for the super- 
hump period by solving the eigenvalue problem of a linear 
one-armed mode in an inviscid disc. Their arguments corre- 
spond with those given above and their results agree very 
well with those of Lubow (1992). They showed that an in- 
crease in disc temperature allowed the eccentric mode to 
propagate further into the disc, with the precession rate be- 
ing reduced as a result. Lubow (1992) also completed several 
two dimensional hydrodynamic disc simulations, in which 
he isolated the various contributing factors to the preces- 
sion. His numerical results clearly showed pressure forces 
to be important. In fact, for those particular simulations, 
the pressure contributions reduced uj to half the dynamical 
value. Till now, observational data has only been compared 
with dynamical estimates for the precession that we know 
to be inadequate. 



3 THE THEORY AND SIMULATIONS 
COMPARED 

The numerical simulations of Whitehurst (1988) lead di- 
rectly to the eccentric disc model for superhumps. In this 
section we compare subsequent numerical results (Hirose & 
Osaki 1990; Kunze, Speith & Riffert 1997; Murray 1998; 
Simpson & Wood 1998; Murray et al. 1999) with the dynam- 
ical equation for precession (equation ^ . We do not attempt 
a comparison with the hydrodynamic equation as the sim- 
ulations of Hirose & Osaki were completed using the sticky 
particle method due to Lin & Pringle (1976) which does not 
account for pressure forces. 

For approximately circular orbits, the rate of dynamical 
precession 



where 



dr 



(5) 
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6s"'' is the standard notation for a Laplace coefficient from ce- 
lestial mechanics. We have made use of the hyper-geometric 
series expression for h\?^ found in Brouwer & Clemence 



1.0 



0.8 




0.10 0.20 



0.30 0.40 
radius 



0.50 0.60 



Figure 1. a{r) is the term encapsulating the radial dependence 
of the dynamical precession rate. Calculated using the hyper- 
geometric series (Brouwer & Clemence, 1961). 



(1961) (equation 42, Chapter 15) to evaluate a{r) as a func- 
tion of radius (see figure hi). Clearly, the precession rate of 
a ballistic particle is an increasing function of radius. But 
as mentioned above, differential precession is not observed. 
Thus in previous applications of the dynamical theory to 
discs, it was assumed that the entire disc precessed as if 
it were a single ballistic particle at the resonance radius, 
Trcs — 0.477 d, (see e.g. Hirose & Osaki). Unfortunately, 
most of these comparisons are of little value because they 
either failed to correctly evaluate a(r) (e.g. Patterson 1998) 
or they made use of an incorrect equation for the precession 
that first appeared in Whitehurst & King (1991). 

Observational data is usually presented in terms of the 
superhump period excess 



Psh 



In terms of the disc precession rate uj, 



(7) 
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So in figure |^ we have plotted as a function of mass ratio 
the superhump period excesses obtained from the disc sim- 
ulations of several different authors. The solid curve shows 
the superhump period excess estimated by equation ^ with 

r — '^rcs. 

All authors used the smooth particle hydrodynamics 
(SPH) technique, except Hirose & Osaki who used the sticky 
particle technique. The pressure contribution to the disc pre- 
cession is thus absent from their simulations. The various im- 
plementations of SPH are described in Flebbe et al. (1994), 
Simpson & Wood (1998), and Murray (1996). 

The calculations completed for Murray (1998), and 
Murray et al. (1999) were of very cool isothermal discs. In 
units of the binary separation, d, and the orbital angular 
frequency, 0"^.^, we set the sound speed c = O.02dr2orb. The 
pressure contribution is proportional to the square of the 
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Figure 2. Collation of superhump period excesses from; Murray 
(1998), Murray et al. (1999) and this paper (filled circles), Hirose 
& Osaki (triangles), Simpson & Wood (asterisks), and Kunze et 
al. (hollow squares). The solid curve is the dynamical contribution 
to the period excess at the resonance radius (as calculated using 
equation ^ from this paper, Simpson & Wood, and Kunze et al. 
lie very close to one another on this plot. 



sound speed and so was very small. We also used a very 
large value for the shear viscosity. This reduced the viscous 
time scale to a value that enabled us to follow the evolution 
of discs to steady state in a reasonable amount of computing 
time. However with a larger shear viscosity it was easier for 
material to penetrate the resonance. Significantly eccentric 
discs and very strong superhump signals were obtained as a 
result. Furthermore, the large shear viscosity inhibited the 
inward propagation of the eccentricity and further reduced 
the effectiveness of the pressure term in equation ^. 

Without the benefit of any significant pressure contri- 
bution to the disc precession, the points from Murray (1998) 
and Hirose & Osaki lie very close to the dynamical precession 
curve. Some are marginally above the curve. As mentioned 
above, the large shear viscosities used allowed the discs to 
effectively penetrate the resonance and extend to somewhat 
larger radii. On the other hand the four points from Mur- 
ray et al. (1999) are all well above the curve, and cannot be 
explained in terms of large discs. For these mass ratios the 
resonance lies beyond the tidal truncation radius. That is, 
for g > 1/4 simply periodic orbits begin intersecting one an- 
other inside rrcs- The orbits in this region will no longer be 
approximately circular, in contradiction of the assumption 
underlying the analytical expressions in Hirose & Osaki, and 
Lubow (1992). We conclude therefore that, gas pressure con- 
siderations aside, expressions such as equation M will only be 
valid for mass ratios less than approximately 1/4. 

The discs of Simpson & Wood, and of Kunze et al. pre- 
cessed significantly more slowly than expected from purely 
dynamical considerations (these points lie below the dynam- 
ical curve in figure ^). Kunze et al. assumed each element 
of the disc radiated as a blackbody and set the sound speed 



according to the blackbody temperature. Simpson & Wood 
used a polytropic equation of state with index 7 = 1.01 
and integrated the internal energy equation for each par- 
ticle. Although a comparison is not straightforward, these 
discs would have been somewhat hotter than the ones con- 
structed for this paper. As a consequence, the retrograde 
contribution to the precession due to pressure forces would 
have been larger, and the values for Pd obtained by Simpson 
& Wood and Kunze et al. are closer to observed P^h- 

We are currently running a set of simulations with the 
shear viscosity ten times smaller than in Murray (1998), but 
with other parameters unchanged. In particular, the sound 
speed c = O.02dr2orb. This gives an effective value of the 
Shakura-Sunyaev parameter a ~ 0.16 at the resonance (as 
opposed to ~ 1.6 in the previous calculations). At present 
we have results only for g = 0.10 (the simulations now take 
much longer to complete). This particular calculation ran 
for 100 orbital periods. As is to be expected, the superhump 
amplitude was reduced. At the conclusion of the calculation 
the superhump period (averaging over 25 superhump cycles) 
Psh = (1.0295 ±0.0005) Porb. This period is significantly be- 
low the dynamical precession curve even though the disc is 
very cool, and it corresponds very well with the results of 
Simpson & Wood, and Kunze et al. Our result has been in- 
cluded in figure ^ but we emphasize that at the conclusion of 
the simulation, the disc mass was still very slowly increasing, 
and the superhump period had not completely stabilised. 

In the above discussion we did not refer to the calcula- 
tions of Whitehurst (1994) simply because he used an initial 
mass transfer burst to set up his discs, and the superhump 
periods were not obviously steady state values. The simula- 
tions described in section 4.2 illustrate the time scale upon 
which resonant discs adjust to changes in the mass transfer 
rate. Before we leave this subject, a word of warning. We 
discovered in one particular trial simulation that once mass 
return from the outer edge of the disc to the secondary star 
occurred, the superhump period excess was reduced by ap- 
proximately 5%. The explanation is simple. Material at the 
outer edge of the disc precesses most rapidly. Remove that 
material and you slow the disc precession. This effect can 
hinder comparison with theory and other simulations. Such 
mass return occurred in the simulations of Hirose & Osaki 
and of Simpson & Wood. 

We conclude that for mass ratios less than 1/4, equa- 
tion ^ provides a useful upper limit for the steady state pre- 
cession rate of a gaseous disc. However, for systems with 
q ^ 1/4, the intersection of simply periodic orbits near the 
resonance radius renders equation n invalid. 

The next step is to compare simulation results with 
eigenmode calculations such as those performed by Hirose & 
Osaki (1993). They tabulated disc precession rates as a func- 
tion of the sound speed at the resonance for systems with 
mass ratios g = 0.15 and 0.05. Both the eigenmode calcula- 
tions and the hydrodynamic simulations show similar mag- 
nitude pressure contributions to the precession. However, as 
yet the assumptions made in making the two sets of calcula- 
tions (e.g. equation of state) differ sufficiently as to prevent 
a more detailed comparison. 
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Figure 3. Observed superhump period excesses (data from Pat- 
terson, 1998) plotted against orbital period. Overlaid are curves 
for the dynamical precession as calculated using equation and 
assuming a white dwarf mass M^^ = 0.76 ± 0.22 Mq, and the 
secondary mass period relationship from Smith & Dhillon (1999), 
i.e. Mscc/Mq = (0.038 ±0.003)P(i-5**±" 0^). The outlying curves 
incorporate the uncertainties in the Smith & Dhillon values. 



4 THEORY AND OBSERVATION COMPARED 

Patterson (1998) tabulated the superhump period excesses 
for 53 systems. In this section we will compare this data with 
both the dynamical and hydrodynamical equations for pre- 
cession. In order to do this we require an independent means 
of estimating a system's mass ratio. Unfortunately, although 
the mass of the secondary star as a function of orbital period 
is reasonably well constrained, there is evidence of consider- 
able variation in the white dwarf masses of otherwise similar 
systems (see e.g. Smith & Dhillon 1998, figure 5). The re- 
sulting uncertainty in q is large enough to interfere with our 
comparisons with theory. An early attempt to compare the- 
ory and observation by Molnar & Kobulnicky (1992), was 
hampered in exactly this fashion (see their figure 2). 

In figure ^ we have plotted the superhump data from 
Patterson (1998) against theoretical predictions for the dy- 
namical precession (made using equation ^ with r = rros). 
In order to estimate the mass ratio of a system of given or- 
bital period we have used the observationally derived stellar 
masses obtained by Smith & Dhillon (1998). They found 
the mean white dwarf mass for all cataclysmic variables to 
be 0.76 ± 0.22 M0, and they estimated that Msoc/Mq = 
(0.038 ± 0.003)P'^ ''*=^° °®^ The central (darker) curve in 
figure ^ represents the best theoretical estimate. The outer 
two curves show the consequences of the uncertainty in the 
Smith & Dhillon values for the theory. In fact. Smith & 
Dhillon found the mean white dwarf mass for systems below 
the period gap to be 0.69±0.13 Mq, and to be 0.80±0.22 Mq 
for systems above the period gap. However precession rates 
recalculated using the adjusted white dwarf masses differed 
only marginally from the original estimates. 

Assumptions about stellar masses aside, figure H shows 



that the superhump observations cannot be adequately ex- 
plained in terms of simple dynamical precession. The retro- 
grade precession due to pressure forces is necessary to bring 
closer agreement between theory and observations. Previ- 
ously published figures showing closer agreement between 
dynamical precession estimates and observations only did 
so because the dynamical precession was calculated incor- 
rectly. For example Patterson (1998) used the equation 



l^dyn/f^orb — 0.23 ■ 



(9) 



The coefficient should be approximately 0.4 (see Hirose & 
Osaki 1990; Lubow 1992; figure ^ above). In order to explain 
superhumps in long period systems with equation ^, mass 
ratios as high as 0.5 are required. These q values are clearly 
incompatible with eccentricity excitation at the 3 : 1 inner 
Lindblad resonance. However, if we use the more defensible 
equation ^ for the precession, then smaller mass ratios are 
obtained. 

Our confidence in the comparison between theory and 
observation is limited by uncertainty in q. The inadequacy of 
the dynamical expression for precession is much more clearly 
apparent when we consider the eclipsing systems OY Car, 
HT Cas and Z Cha. Very accurate determinations of the 
mass ratios of these systems are respectively obtained in 
Wood et al. (1989), Home, Wood & Stiening (1991), and 
Wood et al. (1986). In Table |l| we list for each system the 
observed lu and the dynamical precession rate. In each case 
the difference between the two values is too large to be ex- 
plained by uncertainties in either q or equation ^. The dif- 
ference is simply the (retrograde) pressure contribution to 
the precession. 

We can check whether the inferred cupr is consistent with 
the assumption that the eccentricity is tightly wound. The 
pitch angle i of the spiral wave is given by cot i — \kr\. Thus, 
substituting for k in equation ^ we find that the pressure 
contribution to the precession rate at the 3 : 1 resonance 



--"orb (,7:; T7 -) ■ 
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(10) 



Note that Lubow (1992) had the constant of proportionality 
inverted in his equation 21. If we assume the sound speed 
at the resonance radius c = 0.05 dfloA then we obtain pitch 
angles of 17, 16 and 15° for OY Car, Z Cha and HT Cas 
respectively. These are certainly consistent with the tight 
winding approximation. 

We complete this section with a comparison of the ob- 
servations with the predictions of the hydrodynamic theory. 
As it is not clear what value k should take in equation ^, we 
make the naive assumption that the pressure contribution 
to precession in all 53 systems tabulated in Patterson (1998) 
wiU be similar to that of OY Car, HT Cas and Z Cha. We 
take a mean value for LOpr = — 1.9radday~^, and plot pre- 
dicted precession rates against the observations in figure ^ 
A much improved fit is achieved. As with figure |^, three 
curves are drawn to show the infiuence of uncertainty in our 
knowledge of white dwarf and secondary masses. Much of 
the difference between the curves is due to the uncertainty 
in white dwarf mass, with the uppermost curve correspond- 
ing to Mwd = 0.54 Mq. The observational data is perhaps 
best fit with a curve generated assuming a white dwarf mass 
~ 0.65 M0. 

We recall that equation ^ underestimates the dynamical 
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Table 1. The observed disc precession rate (in rad day~^) for 
three systems with accurately known mass ratios. The dynamical 
precession rate is calculated using equation ^ The pressure con- 
tribution is then simply the difference between those two values. 
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OY Car 


0.102 ±0.003 


1.977 


3.869 


-1.892 


Z Cha 


0.15 ±0.03 


2.961 


4.712 


-1.751 


HT Cas 


0.15 ±0.03 


2.75 


4.773 


-2.02 
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Figure 4. Observed superhump period excesses (data from Pat- 
terson, 1998) plotted against orbital period. Overlaid are curves 
for the hydrodynamic precession as calculated using equation ^ 
with ojpr = — 1.9radday~^. We assume a white dwarf mass 
^wd = 0.76 ± 0.22 M0, and the secondary mass period rela- 
tionship from Smith & Dhillon (1999), i.e. M^cc/Mq = (0.038 ± 
0.003)P'^'^*^'' ''^) . The outlying curves incorporate the uncer- 
tainties in the Smith & Dhillon relation. 

precession rate for systems with mass ratios ^ 1/4. Thus 
the theoretical curves should rise more steeply long-wards 
of say Porb ~ 3 hr. The implication then is that long period 
superhumpers have significantly more massive white dwarfs 
than do their shorter period counterparts. As a consequence 
the mass ratios of the long period systems will be smaller 
than previously estimated. This is qualitatively consistent 
with the eccentricity being excited at the 3 : 1 Lindblad 
resonance, and with our previous result (Murray et al. 1999) 
that the excitation can occur at mass ratios < 1/3. In other 
words, of those systems with say Porb — 3 hrs, only those 
with more massive white dwarfs will exhibit superhumps. 
This then is distinct from any effect caused by systematic 
variation in white dwarf masses in the cataclysmic variable 
population as a whole. 

5 CONCLUSIONS 

We have compared the theoretical predictions for the preces- 
sion rates of eccentric discs with simulation results and with 



observed superhump periods. Comparison with disc simula- 
tions showed that the dynamical equation provided a useful 
upper limit for the disc precession rate of systems with mass 
ratios < 1/4. 

Previous papers showing good agreement between ob- 
servations and a simple dynamical model for precession are 
based on incorrectly calculated dynamical expressions. The 
consequence of using these models is that very large mass 
ratios are predicted for the long period superhumpers. At 
such large mass ratios the resonance lies well beyond the 
truncation radius of the disc and eccentricity cannot be ex- 
cited. 

We show that when the correct expression for dynamical 
precession is used, there is poor agreement with the obser- 
vations. This can be seen even though there is considerable 
uncertainty in estimating the mass ratio of a system of given 
orbital period. 

The observed superhump periods of the eclipsing sys- 
tems OY Car, Z Cha and HT Cas are shown to be signifi- 
cantly less than predicted by the dynamical theory, demon- 
strating that the retrograde contribution of pressure forces 
is important. 

The inclusion of a retrograde pressure contribution to 
the precession rate not only improves the fit to the data but 
also requires long period superhumpers have smaller mass 
ratios than previously thought. This in turn is more con- 
sistent with the eccentricity being generated at the 3 : 1 
Lindblad resonance. 
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